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Abstract 

We examine the behavior of the diffusion coefficient of the ST2 model of water over a broad 
region of the phase diagram via molecular dynamics simulations. The ST2 model has an accessible 
liquid-liquid transition between low-density and high-density phases, making the model an ideal 
candidate to explore the impacts of the liquid- liquid transition on dynamics. We locate charac- 
teristic dynamical loci in the phase diagram and compare them with the previously investigated 
thermodynamic loci. The low-density liquid phase shows a crossover from non-Arrhenius to Arrhe- 
nius behavior, signaling the onset of a crossover from fragile-to-strong behavior. We explain this 
crossover in terms of the asymptotic approach of the low-density liquid to a random tetrahedral 
network, and show that the temperature dependence of the diffusion coefficient over a wide tem- 
perature range can be simply related to the concentration of defects in the network. Our findings 
thus confirm that the low-density phase of ST2 water is a well-defined metastable liquid. 
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I. INTRODUCTION 



It has long been appreciated that water, the most important of all liquids, defies descrip- 
tion in terms of simple liquid behavior in most respects Pfl Many of the anomalies of 
the thermodynamic and transport properties can be attributed to the hydrogen bonds that 
dominate the intermolecular interactions [3-8J. Of the many important studies of water 
conducted over the last several decades, the 1992 proposal that a liquid-liquid phase transi- 
tion (LLPT) occurs in supercooled water has had a particularly significant impact on water 
research [0]. In this proposal, two phases of liquid water, a low density liquid (LDL) and a 
high density liquid (HDL), become distinct below a critical point located in the supercooled 
regime of the phase diagram. H. E. Stanley and co-workers have pursued the implications 
of this proposal over the last 20 years [T0H38] , and much has been learned about the impact 
of a liquid-liquid transition on the properties of water and related systems [391 HO]- 

The LLPT proposal remains controversial because its confirmation via experiments on 
the bulk liquid has been thwarted by rapid ice nucleation at the conditions at which the 
critical point is predicted to occur; bulk studies of the ice melting lines provide indirect 
evidence of a transition [18]. The central strength of the LLPT proposal is that it ratio- 
nalizes the thermodynamical and dynamical anomalies of the supercooled liquid while, at 
the same time, it accounts for the occurrence in experiments of two widely different forms 
of amorphous solid water (low density and high density amorphous ice) as the sub-glass- 
transition manifestations of the LDL and HDL phases [T^ [271 EH I^THH] . Indeed, the 
possibility of a LLPT has now been investigated across the entire class of liquids in which 
tetrahedral bonding dominates the local structure. This class of systems includes water. 
Si [15], Si02 [IS], and BeF2 [H], as well as nanoparticle liquids tailored to exhibit tetrahe- 
dral interactions [HI [H] . Additionally, it has also been shown that liquids with symmetric 
interactions, but with a competition between low density and high density packings, may 
also exhibit LLPT behavior [221 ESI [2H1 ED] • 

Due to the challenges imposed by crystallization on experiments of supercooled water, 
computer simulations have played a central role in the development of the LLPT proposal 
for water, and other tetrahedral liquids. While a LLPT occurs in a variety water models, 
one of the most accessible and clearest examples is in the venerable ST2 model [51j, one 
of the earliest point-charge models for water. "ST2 water" has been extensively studied 
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to clarify the nature of such a LLPT, in particular with regard to thermodynamic and 
structural properties [HI [ISl E2HSZ!- While much is known about the dynamics of ST2 
water P, [Ml EU [SHHS2] a comprehensive study of the dynamical properties comparable in 
scope to the thermodynamic studies is lacking. Therefore, in the present work, we focus on 
the dynamical properties of the ST2 model over a wide range of states that encompasses 
the vicinity of the LLPT. We show that there are striking differences in the nature of the 
dynamics of the HDL and LDL phases. The LDL phase presents a particular challenge, as 
the relaxation time of the liquid increases extremely rapidly with decreasing temperature 
T. However, our results demonstrate that the equilibrium dynamical properties of the 
LDL phase can be understood from the behavior observed in the region accessible to our 
simulations. These results emphasize the central role of the developing network of hydrogen 
bonds for understanding the behavior of the liquid, especially of the LDL phase. 



II. SIMULATIONS 

Our data are generated via molecular dynamics simulations of the ST2 potential [5T] . 
The simulations follow the same protocol as those in Ref. [55] • Specifically, we attempt 
simulations at thermodynamic state points in the density range p = 0.80 g/cm^ tol.20g/cm3 
at intervals Ap = 0.01 g/cm^, and the temperature range T = 200 K to 400 K at intervals 
of AT = 5 K. At each state point considered, we simulate = 1728 molecules, unless 
noted otherwise. We are able to equilibrate the liquid at all chosen densities for T > 255 K; 
for T < 255 K, the lowest density studied is limited by the extreme length of simulation 
required. In total, we examine more than 1500 different state points. The configurations 
used to start simulations are the final configurations from Ref. p3] , where a detailed equation 
of state (EOS) of ST2 water was evaluated. Each simulation is run until the mean-squared 
displacement (MSD) of oxygen atoms reaches 1 nm^ (roughly three interparticle spacings), 
much longer than needed to observe diffusive behavior. We use an integration time step of 
1 fs and employ the Berendsen heat bath with a time constant of 2 ps to control T during the 
production run, in order to compensate for a possible minor drift of the energy in the very 
long simulation runs. The electrostatic potential is truncated at 7.8 A and the energy and 
pressure are corrected using the reaction field method |63] ; periodic boundaries are used to 
minimize size effects. Trajectory information was written to disk every 100 fs for T > 250 K, 
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FIG. 1: (a) Pressure P and (b) potential energy [/ as a function of volume V along isotherms for 
ST2 water, from our N = 1728 simulations. In both panels, we show isotherms from T = 200 to 
350 K, in 5 K steps, from bottom to top. In (a), each isotherm is shifted by cT, with c = 10 MPa/K, 
in order to facilitate comparison of the curves. In (b), the minima of U occurring at low T and 
in the vicinity oiV = 1.2 cm'^/g (p = 0.83 g/cm^) identify conditions at which a particularly well 
developed random tetrahedral network occurs in the liquid. 

and every 10 ps for T < 250 K. 

Figure [T] summarizes the thermodynamic properties of the ST2 liquid in the region ex- 
amined here. Fig. [T]^a) shows the pressure P as a function of volume V along isotherms. 
As T decreases, these isotherms first inflect, and then become progressively flatter, leading 
to the realization of a liquid-liquid critical point at which {dP/dV)T = {d'^P/dV'^)T = 0. 
The critical point conditions in ST2 water have been estimated to occur at Tc = 247 ± 3 K, 
Pc = 185 ± 15 MPa, and pc = 0.955 ± 0.01 g/cm^ |56l M\- 

Figure [l](b) presents the variation of the potential energy [/ as a function of V along 
isotherms of the liquid. As noted in previous work |S31 El] , the emergence with decreasing 
T of significant negative curvature in these isotherms is a thermodynamic signature of the 



4 



approaching liquid-liquid instability: To be stable with respect to phase separation, the 
Helmholtz free energy A must have positive curvature {d'^A/dV'^)T > 0. Since A = U — TS 
(where S is entropy), the negative curvature shown in Fig. [T]^b) indicates an energetic driving 
force for phase separation that becomes dominant for T < Tc [53j. It is also interesting 
to observe the minima of U occurring at low T and in the vicinity of ^ = 1.2 cm^/g 
(p = 0.83 g/cm^). The presence of this minimum has also been observed in several other 
models of tetrahedral liquids [53 | IMI - 166] and signals the optimal network volume (or density), 
i.e. the conditions at which a particularly well developed random tetrahedral network occurs 
in the liquid. Indeed, hydrogen bonding requires a well-defined distance and orientation to 
be effective, and at the optimal density (which is usually close to the fully bonded open 
crystal density) geometric constraints allow the establishment of a fully bonded disordered 
network. 



III. DIFFUSIVITY OF THE ST2 MODEL 

For each state point simulated, we calculate the mean-squared displacement (r^(t)) and 
evaluate the diffusion coefficient from the Einstein relation 

D = hm (1) 

Here (...) represents an average over possible time origins. Since we have more than 1500 
state points to consider, we automate the evaluation of D by a linear fit of (r^(t)) for all data 
such that (r^(t)) > 0.5 nm^, a restriction that ensures that all fitted data are well within 
the diffusive regime. 

We first consider D along isotherms, as shown in Fig. [2j For T < 335 K, D exhibits a 
(weak) maximum with increasing p or P, as known experimentally [671 168]. This feature is 
normally attributed to the breaking of hydrogen bonds with increasing density, which allows 
for increased diffusion, until packing considerations become dominant and D decreases. As 
compared with water, the ST2 model overestimates the pressure of the maximum in D, 
which is not surprising given that ST2 overemphasizes the tetrahedral structure relative to 
water [69]. Hence, there exists a locus of points (which we denote -Dmax) on the EOS surface 
at which D is a maximum along isotherms; the shape of this locus in the phase diagram 
is discussed below. Fig. [2] also includes data obtained for T < T^. The filled data points 
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FIG. 2: The diffusion coefficient D along isotherms as a function of density p ((a) and (b)) or 
pressure P ((c) and (d)). The lines are intended only as a guide for the eye. We plot only a 
selection of our available isotherms to avoid overcrowding in the figure. We show both log ((a) and 
(c)) and linear ((b) and (d)) scales to cover a broad range of D (log scale) and to show the form 
of the maxima of D (linear scale). The solid symbols indicate state points in the unstable zone of 
the liquid-liquid coexistence region. Note that L> is a multi-valued function of P along isotherms 
for T < 250 K because P is a non-monotonic function of p in the co-existence region along these 
isotherms. 

indicate simulations in the unstable regime where our system (simulated at constant volume) 
is phase separated into regions of LDL and HDL. The values of D in this region thus reflect 
a weighted average over the LDL and HDL phases. We note that D decreases by nearly two 
orders of magnitude as the system progressively transforms from pure HDL to pure LDL 
along the lowest T isotherms. 

We next present the T-dependence of D along isochores and isobars in Fig. |3| At high 
T, D is described by the expected Arrhenius behavior 

D = Doexp[-E/kBT], (2) 
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FIG. 3: Diffusion coefficient D as a function of T along (a) isocliores and (b) isobars. Along 
isochores the solid symbols indicate systems from the unstable zone of the coexistence region. No 
data for unstable states are shown along the isobars in (b), since D is a multi- valued function of 
pressure in this region [see Fig. [2]j^c,d)]. We show only a selection of densities and pressures to 
avoid overcrowding the figure. 

where E is the activation energy for diffusion and Dq is the hmiting high-T diffusion coeffi- 
cient (both determined from fitting eq. |2]to the data). On coohng below a temperature Ta, 
D exhibits so-called "super-Arrhenius" behavior, where D decreases faster than expected 
relative to the high-T behavior. This rapid decrease is typical of glass-forming liquids as 
they approach the glass transition temperature Tg. We estimate Ta for D along either iso- 
chores or isobars by finding that T at which ksT / E \n{D / Dq) > 1.02; by construction, this 
quantity must equal unity for high-T Arrhenius behavior |70] . 
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FIG. 4: The difFusivity exponent 7 as a function of (a) density and (b) pressure. The pressure 
dependence of 7 parallels that of water in the region where experimental data is available. Exper- 
imental data are from Ref . \<o7\ [68] . The lines are only a guide to the eye. 

The non-Arrhenius behavior of D for many glass-forming liquids is well-accounted for by 
the mode-coupling theory (MCT) for the glass transition [71], which predicts that 

D ~ (T - Tmct)^ (3) 

where Tmct is the temperature of an avoided vitrification of the ideal MCT. Although discus- 
sion continues concerning the region of validity of the MCT approach, the phenomenological 
appearance of power-law behavior immediately below Ta is widely acknowledged. Moreover, 
power-law behavior of D in water at atmospheric pressure has been appreciated since the 
1970's [72]. This power-law behavior of D in water was first interpreted as possible evidence 
for an underlying spinodal singularity, but more recently it has been connected with the 
MCT approach [731 - I76] . Accordingly, we fit separately the isobaric and isochoric data for D 
to eq. [3] to evaluate the locus of Tmct in the phase diagram, which we discuss below. Fits 
must be considered along both isochoric and isobaric paths, since the diffusivity exponent 
7 will depend on the path of the approach to Tmct- Additionally, care must be taken when 
making this fit, since it is apparent from Fig. [3] that it appears D may return to Arrhenius 
behavior at lower T, a phenomenon which we will discuss in detail in the following section. 
Therefore, we have excluded data where the behavior may revert to Arrhenius. 

We find that the diffusivity exponent 7 is non-monotonic with density, reaching a max- 
imum of ~ 2.8 at p ~ 0.88 g/cm^ (Fig. |4]). The exponent 7 is also non-monotonic as a 
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FIG. 5: The (a) T — p and (b) P — T phase diagrams, combined with the dynamical features 

('X') 

determined here: -Dmax! the locus of extrema of D along isotherms; , the onset of non-Arrhenius 
behavior; T^, the breakdown of Stokes-Einstein behavior (from [62l|78]); TmcT' *he extrapolated 
mode-coupling divergence temperature. The superscript X for the last three loci indicates whether 
the locus was determined along isochores {X = p) or isobars {X = P). The color map indicates 
the value of logD; red (blue) represents the highest (lowest) values of D. The liquid-gas spinodal is 
indicated by red diamonds; the liquid-liquid spinodals are indicated by blue triangles. The liquid- 
liquid critical point location is indicated by the open circle. The locus of specific heat maxima 
(jiaax ^j-g^gjjg ^j^g maxima of Cp along isobars; the dotted portion of the line is an extrapolation to 
the critical point, where Cp diverges, pext labels the locus of density extrema; inside the pext locus, 
the isobaric expansivity is negative, while outside it is positive. All thermodynamic data are from 
Ref. [55]. 

function of pressure. A comparison with experimental data for 7 shows that the behavior of 
7 for ST2 is roughly parallel in the region where both experimental and simulation data are 
available. This is in contrast with the SPC/E model of water, where the opposite pressure 
dependence occurs [76]. 

To summarize the properties of we collect the resulting characteristic features from 
the isothermal, isochoric, and isobaric plots (Fig. |2]and Fig. [s]) and plot them together with 
the known thermodynamic features. We show these features, along with a colormap for D, 
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in the p — T plane and in the P ~ T plane in Fig. |5} As expected from Ref . [77] , the locus 
-Dmax hes outside the region of negative entropy-volume correlations, delineated by the locus 
of extrema of density pext- The existence of a maximum in D also results in non-monotonic 
behavior of Ta and Tmct, since these temperatures represent a nearly constant value of D. 
The locus of Tmct is nearly coincident with the lower bound of our simulated data, since this 
also represents the time scale where simulations become prohibitively lengthy. We also show 
the locus Tx of the breakdown of the Stokes-Einstein (SE) relationship from Ref. |62t 178] . It 
is interesting to note that Tx > Tmct, since the breakdown of the SE relation is often linked 
with Tmct- For low pressure, Tx is correlated with (but not coincident with) the extrema 
of the specific heat [78]. Finally, the shape of the loci Ta, Tmct and Tx are all similar, with 
the non-monotonic behavior becoming more pronounced with lower T. 

IV. STRUCTURE AND DYNAMICS OF THE LOW DENSITY LIQUID PHASE 

One of the most challenging aspects of the liquid-liquid phase transition in ST2 water 
is characterizing the properties of the LDL phase. As shown in the previous section, at 
low density {p < pc or P < P^.) the value of D decreases far more rapidly as T decreases 
than for p > p^. Obtaining equilibrium values for both structural and dynamical properties 
of the LDL phase is therefore particularly demanding of computational resources. In this 
section, we examine the behavior of the liquid in this low density regime, and show that 
many properties of the LDL phase can be determined from the behavior observed in the 
region accessible to our simulations. 

In particular, we focus on states along the optimal density p = 0.83 g/cm^ isochore. 
As shown in Fig. [l|^b), this density approximately corresponds to a minimum in isotherms 
of U versus V, indicating that the structure of the random tetrahedral network (RTN) 
is particularly well developed at this density. We therefore expect that the characteristic 
properties of the LDL phase will be most prominent at this density. 

To complement the = 1728 simulation results described in the previous section, we 
examine an ensemble of 40 independently initialized and equilibrated simulations along the 
p = 0.83 g/cm^ isochore for a system of size A^ = 216. The smaller system size allows us to 
probe longer time scales than for A^ = 1728. In addition, by averaging our results at each T 
over the 40 independent runs, we can significantly reduce the statistical error of our results. 
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FIG. 6: Arrhenius plots of D, a • fj, and b ■ (eis — Cjg)^; fd is the fraction of defects found in 
instantaneous liquid configurations. All curves are for liquid ST2 water along the p = 0.83 g/cm^ 
isochore, as determined from our N = 216 simulations. D is plotted in cm^/s, while /J and 
(eis — Cjg) are multiplied by arbitrary constants a and b to facilitate comparison within the same 
plot. The straight lines are fits to an expression proportional to exp{—Ea/RT) for the four lowest 
T points along each curve; each of these lines is labeled by the activation energy Ea obtained from 
the fit. 

We thus use these smaller systems to carefully parse the behavior of the low-density liquid. 

To evaluate the dynamical behavior of the = 216 system at p = 0.83 g/cm^, in each of 
the 40 runs, the diffusion coefficient is estimated from D = (r^)/6t, where (r^) is the mean 
squared displacement at the end of the run, and t is the time of the run. All our production 
runs for = 216 are carried out until (r^) = 1.0 nm^, or t = 0.5 ns, which ever takes longer 
to achieve. At the lowest T (255 K), the longest runs require up to 350 ns. We then average 
the results over the 40 runs. 

The results for D as a function of T are shown in an Arrhenius plot in Fig. [6j Consistent 
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FIG. 7: Behavior of liquid ST2 water as a function of T along the p = 0.83 g/cm^ isochore, 
determined from the = 216 simulations, (a) U — 3RT for the liquid (circles), and for crystalline 
ice I/i (solid line), compared with the behavior of eis (squares) from the inherent structures. The 
dotted green line shows the estimated limiting ejg for amorphous states, larger than the ice energy, 
(b) The isochoric heat capacity Cy of the liquid. The solid line is only a guide for the eye. (c) 
The fraction of defects fd as found from instantaneous (circles) and inherent structure (squares) 
configurations of the liquid. 

with the results of Fig. |3} we find that the T-dependence of D crosses over from non- 
Arrhenius (fragile) at higher T to Arrhenius at the lowest T with an activation energy 
E = llA kJ/mol. This appears to be the beginning of a "fragile-to-strong crossover" (FSC); 
such behavior has been observed and studied in a number of systems in which RTN-like 
structure emerges at low T [25|, EU |65l [TSj ITOl - IST] . Since we have not yet reached the 
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low energy RTN, the activation energy has not reached its low-T asymptotic value. The 
approach to the RTN can exhibit intermediate Arrhenius behavior over several decades in 
D before reaching the asymptotic limit [25] . 

To support our interpretation that we are approaching strong liquid behavior as T de- 
creases, we test the expectation that the FSC is associated with the approach to the lowest 
lying minima of the liquid's potential energy landscape (PEL). To confirm this behavior in 
ST2 water, we carry out conjugate gradient quenches of a large number of configurations 
(at least 400 for each value of T studied) in order to evaluate the average energy of a liquid 
configuration when quenched to its nearest local minimum of the PEL, referred to as the 
average inherent structure energy, eis. 

In Fig. [7](a) we show both U and eis as a function of T for p = 0.83 g/cm^. In both 
cases we find that an inflection occurs in the vicinity of T = 285 K, and that at lower T the 
rate of decrease of the energy slows with decreasing T. Fig. [7]^a) also shows U for the ice Ih 
crystal, the value of which sets a lower bound on ejs for the liquid. The relative shape and 
position of these curves suggests that the energy of the low density liquid is approaching a 
low T limit associated with the "bottom" of the liquid PEL that lies above the crystal state. 

The approach of the liquid to the bottom of the PEL is also confirmed in Fig. |8j which 
plots the sampled probability distribution of the set of eis values obtained at low T. While 
at high T the distributions are approximately Gaussian, at the lowest T = 255 K the 
distribution has become distinctly skewed and narrower, reflecting the approach to a finite 
lower bound for the possible values of eis- We assume that the approach of eis toward its 
minimum value ejg obeys 

eis = e?s + i?e-^-/^^ (4) 

and estimate ejg by fitting eis to this form over the range of T in which D appears to follow 
an Arrhenius behavior. We approximate this range as the lowest four T shown in Fig. |6| 
and obtain e^g = —57.6 kJ/mol. As expected for a liquid, this value is above that of the 
value of U for the crystal, ice Ih- The limiting strong behavior will only be apparent when 
the system nears e^g. We notice that the value of ejg is almost exactly half of the value of 
the activation energy for diffusion, supporting the view that D is strongly related to the 
single particle energy. Since each bond connects two molecules, we can evaluate the binding 
energy per bond Eb ~ ejg/2 = —28.8 kJ/mol, which is larger than the bare HB energy 
Ehb ~ 20 kJ/mol [58l|59], since Eh includes non-bonded interactions. 
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FIG. 8: Probability distributions of eis for the liquid from the = 216 simulations, along the 
p = 0.83 g/cBo? isochore. From left to right, T = 255 to 275 K, in 5 K steps. 

We also note in Fig. [7](b) that a maximum in the isochoric heat capacity Cy occurs in the 
vicinity of T = 285 K, and as expected, coincides with the inflections of U and eis shown 
in Fig. [7|^a). This thermodynamic signature emphasizes that even though all our data at 
p = 0.83 g/cm^ lie above T^, for T < 285 K we have crossed into the regime in which the 
properties of the liquid are increasingly dominated by the LDL phase, which is distinct from 
the HDL phase for T < Tc. 

In order to provide some structural insight into the dynamical behavior along the p = 
0.83 g/cm^ isochore, we examine the role of defects in the RTN structure of the liquid. 
In a perfect RTN, all molecules would have exactly four nearest neighbors (nn's) within a 
distance given by the first minimum of the oxygen-oxygen radial distribution function. This 
distance is approximately 0.35 nm in ST2 water at p = 0.83 g/cm^ [SS]- We thus define 
the fraction of RTN defects fd as the average fraction of molecules in the system that do 
not have four nn's within a distance of 0.35 nm; hence such defect molecules will either 
have more than, or less than, four nn's. fd is plotted as a function of T in Fig. ^c), as 
obtained from both instantaneous liquid configurations, as well as from inherent structures. 
In both cases we find that fd decreases as T decreases, passes through an inflection near 
T = 285 K, and approaches zero at our lowest T. Indeed, as shown in Fig. [6} we find that 
the approach of fd to zero at low T closely follows an exponential decay, and in particular, 
that the defect activation energy Ea obtained by plotting is approximately equal to the 
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FIG. 9: Parametric plot of D versus fd (squares) for the liquid from the N = 216 simulations, 
along the p = 0.83 g/cm^ isochore, from T = 255K (bottom left) to 360 K (top right). In this 
plot, fd is evaluated from the instantaneous configurations of the liquid. For each T, we also show 
the "cloud" of points giving the values of D and fd obtained from each of the 40 independent 
simulations conducted at each T. The straight line has a slope of 4. 

value of Ea obtained by fitting D to an Arrhenius law at low T. Since Ea ~ 4i?b, fd is 
approximately proportional to the exponential of Eb/RT; hence, fd essentially measures the 
probability that a single bond is broken. 

To test if the relation between D and fd holds at higher T (where the behavior is non- 
Arrhenius) we present in Fig. [9] a parametric plot of D versus fd over the entire T range 
studied a.t p = 0.83 g/cm^, from T = 360 K to 255 K. Except for some deviations at the 
highest T, we find a remarkably consistent behavior in which D varies as a power law in 
fd over almost four decades, and with an exponent very close to 4. This result strongly 
suggests that the hquid at this density can be understood as a disrupted RTN, and that the 
localized excitations of this RTN {i.e. the defects) control the transport properties of the 
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liquid. 

The relation between D and fd can be anticipated by works on colloids and nanoparticles 
with highly directional interactions. In particular, it was found that in systems of patchy 
colloids with four sticky-spots in a tetrahedral geometry [M] and in nanoparticles linked by 
DNA with tetrahedral orientation [821 ES] D is given by the fourth power of the fraction of 
broken bonds, D ~ One possible explanation for this behavior is that D is controlled by 
the diffusion process of the particles which have four broken bonds, and which are thus free 
to wander around in the available empty space, searching for the rare free dangling ends — 
i.e. for the rare sites of the network with incomplete bonding — to stick. Indeed, if this is 
the case, D is proportional in a first approximation to the fraction of un-bonded particles, 
and this fraction scales with It is important to remember that liquid water at this 
density belongs[83] to the category of so-called empty liquids|83, i.e. liquids in which the 
fraction of space occupied by a space-filling representation of the molecules is significantly 
smaller than the close packing value. For the case of liquid water, the corresponding effective 
packing fraction, when molecules are considered as hard-spheres of diameter 0.28 nm (the 
mean hydrogen bond length), is 38%. 

The comparison between water and these tailored tetrahedral systems [Ml l82| 183] can 
also help clarify the connection between the activation energy of the dynamics and the 
bonding energy. Indeed, in these systems, it is clear that the activation energy for diffusion 
is given by four times the bond energy, since the interactions are short-ranged and there 
is no energy contribution from non-bonded neighbors. For ST2, each hydrogen bond has 
a strength of -Ehb ~ 20 kJ/mol [581 EH], significantly less (by 9 kJ/mol) than the overall 
binding energy per bond Ei, at the lowest T simulated. If breaking of HBs are the limiting 
factor in diffusion, we can anticipate that the asymptotic activation energy for diffusion 
may decrease to a limiting value of ~ 80 kJ/mol, similar to the low-T value (74 kJ/mol) 
estimated for water [25] . 

V. SUMMARY 

In conclusion, we have provided a comprehensive survey of the diffusive properties of the 
ST2 water model in the vicinity of the liquid-liquid transition. Our results demonstrate 
that the structural and dynamical properties of the LDL phase that becomes a thermody- 
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namically distinct liquid phase for T < Tc are already well established in the liquid at low 
density for T > Tc. The LDL phase is here revealed as a highly structured liquid, whose 
properties are dominated by the progressive emergence of a RTN structure as T decreases. 
In this sense, the properties of the LDL phase are entirely consistent with those of low 
density amorphous ice, as determined in experiments. Finally, we show that the dynamics 
of the LDL appears to be fully controlled by the presence of defects of the network, whose 
concentration is controlled by the bond energy. 
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